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5^ ■ Extended calculations of the deuteron's static properties, based on the numerical solution of the 

Bethe-Salpeter equation, are presented. A formalism is developed, which provides a comparative 
analysis of the covariant amplitudes in various representations and nonrelativistic wave functions. 
The magnetic and quadrupole moments of the deuteron are calculated in the Bethe-Salpeter for- 
malism and the role of relativistic corrections is discussed. 

I. INTRODUCTION 

A theory applicable for studying nuclear phenomena, involving momentum transfers of a few GeV or higher, must be 
relativistic. A traditional approach to processes with nuclei, based on the nonrelativistic Schrodinger wave functions, 
is not adequate if a large momentum transfer "sneaks" into the nuclear amplitudes, and the corresponding nucleon 
momentum p becomes large, say p > m (m is the nucleon mass). One can extend the usage of nonrelativistic wave 
functions by incorporating successively the relativistic corrections ~ (p/m) n , however, it will eventually fail at some 
point of p. On the other hand, the nonrelativistic approach was the only one which allowed for the detailed description 
of the static properties of the nuclei and low and intermediate energy nuclear reactions. 

In recent decades, extensive studies of few- nucleon systems were performed within Lorentz invariant models [Q-^). 
The success of these elaborate studies allows one to conclude that the covariant approach has now the capability 
to replace, at least for few-nucleon systems, the approaches relying on nonrelativistic wave functions |(|. Most of 
the phenomenological success in the relativistic treatment of few-nucleon systems has been achieved within such 
models which are based on the covariant meson-nucleon theory and dynamical equations In these models, the 

satisfactory results have been obtained for the nucleon-nucleon scattering, the properties of the lightest nuclei, various 
electromagnetic and hadronic interactions with nuclei, and some advance have been achieved for many-body nuclear 
systems (see e.g. discussions and references in ref. ||). 

The deuteron, as the simplest nuclear system, is an appealing object to be described by the models invented in the 
realm of nuclear physics. At the same time, there is a fair amount of experimental information available about the 
deuteron's properties themselves and reactions with the deuteron. More interesting and precise data is expected after 
the start of the exciting research program at CEBAF. Therefore, there is a possibility to compare exact theoretical 
results with the experimental data in a clear way, not dimmed by extra effects, such as the "more-then-two-body" 
calculations. 

Still, the relativistic approach to the deuteron is not as popular as the one utilizing nonrelativistic wave functions 
|^,^| . There are seemingly two main reasons for this. First, the deuteron, as any other nucleus, is essentially a 
nonrelativistic system, since it is composed of weakly bound massive nucleons. The bulk of the static properties of 
such a system obviously can be fitted in the nonrelativistic approach by adjusting the phenomenological potential or 
the wave function. Besides, the experimental data for the reactions with the deuteron is also mainly available in the 
nonrelativistic domain. Second, the relativistic models, especially those based on field theory, are technically more 
difficult and have a more sophisticated physical interpretation than the nonrelativistic approach. Both these reasons, 
together, define the typical pattern for the attempts devoted to promote the relativistic description of the nuclei. 
These works are usually highly specified for the particular reactions or kinematic domains where the advantage of the 
covariant approach can be explicitly displayed. They are often filled with technical details uncommon for that part 
of the scientific audience which is not directly involved in this research direction. That is why this is so important to 
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have simple and intuitively clear interpretations of the relativistic calculations, and an explicit systematic method to 
compare the relativistic and nonrelativistic results. 

In the present work we are going to analyze the extended calculations of the static properties of the deuteron 
utilizing the Bethe-Salpeter (BS) amplitudes which are recently computed numerically |9|. The main goal of our 
paper is to contribute to the development of the physical intuition for understanding the relativistic calculations and 
their comparison to the nonrelativistic calculations. Our basic idea is to compute the observable densities of various 
charges (e.g., vector and axial- vector charges) in both the relativistic and the nonrelativistic formalisms and use these 
densities as tools to compare relativistic amplitudes and nonrelativistic wave functions, which can not be rigorously 
interrelated otherwise. In doing so we pursue, in some sense, the goals opposite to the ones we outlined above as 
typical for the approach within the covariant description of the deuteron. Another goal of our paper is to fill some 
gap in the literature by giving explicit expressions relating the BS amplitudes in different representations, which will 
help to compare the relativistic amplitudes computed in different models. 

We calculate here the magnetic and quadrupole moments of the deuteron within the Bethe-Salpeter formalism. 
The investigation of these static characteristics of the deuteron is still an important topic in nuclear physics. In 
the nonrelativistic models it gives the direct information about the tensor components in the NN interaction and the 
magnitude of the D wave probability in the deuteron. However, there is an essential problem in fitting the experimental 
values of the quadrupole and magnetic moments within the same D wave probability in the nonrelativistic calculations 
(cf. ref. ]To| and references therein). The efforts, aiming to solve this difficulty, go in two main directions, namely 
calculating the corrections of the meson exchange currents and taking into account the relativistic effects 

p|,p| Jl^Jl5| ,p|, ^6| , ^8| . In the conventional approach the mesonic degrees of freedom and relativistic effects are treated 
as corrections to the nonrelativistic potential theory. It is found that, by adding these effects to the quadrupole 
moment, a satisfactory description of the data may be achieved for a broad range of different potentials (T^| , while the 
magnetic moment shows an essential sensitivity to the model calculations of the meson exchange currents. Moreover, 
the consistency of such calculations is not at all clear. For this reason a comprehensive covariant investigation has 
its own right. A prominent feature of the relativistic consideration within the Bethe-Salpeter formalism is that the 
meson exchange effects due to pair creation currents is taken into account consistently |l0| , |l7| , ^9| , so that the essential 
part of the mentioned effects may be estimated in a self consistent way. 

The general approach to calculate the static characteristic of the deuteron within the BS formalism has been 
elaborated by several authors since some time (see, for instance, refs. P,p^- ^6pl| ] ) and numerical estimates have been 
performed. However, explicit calculations have been done within additional approximations for the solution of the BS 



equation, that is, for a separable interaction and by disregarding the negative energy states |15|, or with one nucleon 
on mass-shell j|, or from a general point of view with adjusting the probability of the P states as to fit simultaneously 
both the quadrupole and magnetic moments ]l6| ] (for this goal one needs an anomalously large pseudoprobability of 
the P waves, say ~ 1.5%). In this paper we perform a covariant calculation of the quadrupole and magnetic moments 
of the deuteron within the exact solution of the BS equaion and avoid additional approximations to the problem. 

Our present investigation is also partially motivated by the renewed interest in the experimental investigation 
of the nucleon and deuteron spin-dependent structure functions at low momentum transfer Q 2 |22| . This interest 
is connected to the study of the Q 2 -evolution of the Gerasimov-Drell-Hearn sum rule piij] , which relates the spin- 
dependent structure functions of the targets to their the magnetic moments. For instance, only a correct description 
of the deuteron magnetic and quadrupole moments will assure a reliable extraction of the information about the 
neutron structure function from the deuteron data. 

Our paper is organized as follows. In Section II the basic covariant formulae for the electromagnetic current and 
static moments of the deuteron are presented. In Section III the general definitions of the Bethe-Salpeter amplitudes 
for the deuteron are given in different representations and their symmetry properties are studied in detail. The 
transformation matrix relating the amplitudes in different representations is determined. The relativistic amplitudes 
are compared to the nonrelativistic wave functions, using the calculated observables, e.g., the vector and axial charge 
densities. In Section IV the covariant formulae for the magnetic and quadrupole moments are derived in the Breit 
frame. The effects of the Lorcntz deformation and dependence of the amplitude on the relative energy of the two 
nucleons in the deuteron are explicitly taken into account. The terms corresponding to the nonrelativistic expressions 
for the moments are determined in explicit form and the relativistic corrections are computed. The Sections V and 
VI contain conclusions and the summary, respectively. 

II. RELATIVISTIC KINEMATICS OF THE ELECTROMAGNETIC CURRENT 

The definition of the quadrupole moment Q r> and the magnetic moment fip of the deuteron appears most trans- 
parent if one starts with the famous Rosenbluth formula for the elastic scattering of electrons off the deuteron, 
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A(q 2 ) + B(q 2 )t&Tl 2 - 

with the following decomposition of the electromagnetic form factors 

A{ q 2 ) = F 2 { q 2 ) + lv 2 F 2 Q { q 2 ) + \vFUe 
B(q 2 ) = ^(l + V )Fi 4 (q 2 ), 



(1) 

(2) 
(3) 



where r) = —q 2 /4M E and M E is the deuteron mass. Q 2 = —q 2 denotes the momentum transfer. Then the quadrupole 
and magnetic moments of the deuteron arc defined via the normalization conditions for the charge (Fc), quadrupole 
(Fq) and magnetic (Fm) form factors at vanishing momentum transfer q 2 = 



F C (0) = 1, f q (o) = m z d q d , f m (o) = m 



M D 



(4) 



with m as nucleon mass. 

The general form of the deuteron electromagnetic current which is invariant under the Lorentz and time-reverse 
transformations, is given by 



(P',A'|JJP,A) 



2Mr 



;(P',A')Jr^(P,A), 



(5) 



where £*(P', A') and e(P, A) are the polarization four vectors of the initial and final deuteron states. The covariant 
normalization of the current reads 



lim (P',A'|J M |P,A) = < 



P 



M S X , X . 



q 2 ^0' M D 

The matrix element can be expanded in terms of the scalar form factors in the form 



3%, = (P' + PY [gpMq 2 ) - ^F 2 (q 2 ) j + {g*q a - g^G^q 2 ). 
The scalar form factors Pi. 2 and G\ are related to the form factors Fc.q.m, by (cf. p3) 

F c (q 2 ) = F^q 2 ) + ^[F^q 2 ) + (1 + V )F 2 (q 2 ) - G^q 2 )}, 

F Q (q 2 ) = F^q 2 ) + (1 + V )F 2 (q 2 ) - G^q 2 ), 
F M (q 2 ) = G^q 2 ). 

In the nonrelativistic impulse approximation these deuteron form factors read 

F C (q 2 ) = (G E (q 2 ) + G E (q 2 ))C E (q 2 ), 
F Q (q 2 ) = (G p E (q 2 )+G*k(q 2 ))C Q (q 2 ), 

F M (q 2 ) = ^ 



(G E (q 2 ) + G n E {q 2 )) C S {q Z ) + ^ (G p M tf) + G n M {q')) C L {q*) 



(6) 



(7) 



(8) 
(9) 
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(11) 
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where G^ E n \q 2 ) ^G^'™' ) (g 2 )^ are the electric (magnetic) nucleon formfactors and the invariant functions C(q 2 ) are 
defined by 



C E (q 2 )= f(u 2 + w 2 )j (^)dr, C E (Q) 


3V2 



= 1, 



2r] 



2^2 
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dr, C Q (0)=M D .Q D , 



(14) 



(15) 
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Cl(<1 2 ) = ljw 2 [jo (f ) +h (f )] dr, G L (0) = |pb, (16) 



oo oo 

c s(a .,-/(^-^)*(?)* + /(^ + ^) *(?)*, 

o o 

G s (0) = 1-|p d . (17) 

Here j, is the modified Bessel function of i-th order, u and w represent the S and D waves of the nonrelativistic 
deuteron wave function, and Pjj is the weight of the D wave in the deuteron wave function. 

To calculate the form factors Fc,q,m within the Bethe-Salpeter formalism one has to express the current (|J) in 
terms of the BS amplitudes and, then, to extract the coefficients of different Lorentz structures given by eq. (Q). 
Taking the limit q 2 — > 0, the static moments can be also obtained. Apparently, these calculations can be done in any 
particular reference frame. For example, the Breit frame is especially convenient for such type of calculations. The 
Breit frame for the deuteron is defined by 

P = -~, P' = f, Po = P^=E. (18) 

Choosing q along the positive z axis and contracting J^ a and the polarization vectors, which obey e(P' , A'), e(P, A): 

e(P',l)= £ (P,l) = --i=(0,l,z,0), (19) 

e(P', -1) = e(P, -1) = -J=(0, 1, -i, 0), (20) 



£ (p,o) = (-v^,o,o,yiT^), £ (P',o) = (V?/,o,o,yTT^), (21) 

one arrives at the expressions for the matrix elements of the deuteron electromagnetic current in terms of the form 
factors: 



(P', A'| J°\P, A) = e^/T+ri (F 1 5 XXI + 2 V [Fi + (1 + r])F 2 - GJ (Wa,o) , (22) 
(P', A'| J X \P, A) = e ^yTT^Gi («v,a+i - Sx>,\-i) , (23) 

(P', X'\J y \P, A) = -ie^-^T+^Gi (<5a',a+i + 5v,A-i) , (24) 

(P',A'|J Z |P,A) = 0. (25) 
Thus the magnetic and quadrupole form factors of the deuteron are recovered by 

m /7T,. (P',A' = l|J a: |P,A = 0) 

^ 2^71+^ V ' 

Equations (^)-(^) are the basic relations providing the calculations of the electromagnetic characteristics of the 
deuteron. In practice, one needs to define explicitly the operator of the electromagnetic current and calculate its 
matrix elements with the deuteron states |P, A). 

III. THE BOUND STATE WAVE FUNCTION 
A. General definitions 

Using the technique presented in ref. j^], the BS equation for a bound state in ladder approximation can be written 
in the form 



4 



tffo,, p) P) + ^ — J dp (p _ p , )2 _ M| = 0, (28) 

K( Po , p) = (El -v\- Ml/ A) 2 - plMl, (29) 

where x(p;P) is the BS amplitude for the deuteron in the matrix representation A(pi) — pi ~ m; p = (po,p) is 
four-momentum of the i-th nucleon expressed in terms of relative four-momenta p or p' and the center-of-mass (cm.) 
momentum P = (A/d,0): pi 2 = P/2 ± p; B enumerates the exchanged mesons 7r, er, w, p, 5, 77; is the mass of the 
meson; is the interaction vertex between the nucleon and the corresponding boson B and Xb = g\ / '47r with 
being the coupling constant. We use here the short hand notation p = p M 7^. 

Since the BS amplitude x an< 4 its adjoint x satisfy the homogeneous BS equation they are determined up to an 
arbitrary constant which is fixed by additional normalization condition. In the ladder approximation the normalization 
constant may be fixed by computing the matrix element of the electromagnetic current at q 2 = 

Tr{x(p;P)^ X (p;P)(m-Pi)} = 2P»- (30) 



The normalization condition (|30| ) coincides with the one used in ref. |14f| . 

The BS amplitude is a (4 x 4) matrix in the spinor space, and consequently the BS equation ( p8| ) possesses this 
matrix structure as well. To solve this matrix equations one can utilize a decomposition of the BS amplitude over a 
complete set of (4 x 4) matrices and solves a system of coupled equations for the coefficients of such a decomposition. 
The choice of the representation of the matrices depends on the concrete attacked problem. Certainly, different 
representations are related by linear transformations, and it is straightforward (but cumbersome) to transform results 
from one representation to another one. In our opinion, to solve the BS equation and to compute matrix elements of 
the deuteron observables (as for instance eq. (^0)) a convenient way is to decompose the amplitude in terms of the 
complete set of Dirac matrices, which form the Clifford algebra (for more details cf. ref. ||). By exploiting the parity 
invariance of the BS amplitude 



Pxd(po,p) = VvloXD(po,-pho, (31) 
it may be written for the deuteron, which has positive parity eigenvalues {rfp = 1), as 

XD (p; P) = 75 P + 7 5 7°A° - (7- V) - 75 (7 -A) - 2*7°( 7 .f °) - 2 7 7 5 (7-T), (32) 

with pseudoscalar (P, A ), axial (A) and vector (T°, T, V) functions depending only upon the relative four momentum 
p in the cm. frame. The angular dependence of the state with spin J = 1 and its projection M. owing to the rotational 
invariance of eq. ( p8| ) is expressed in terms of the spherical and vector spherical harmonics. For example, when denoting 
X = (P, A ) and X = (A,f°,T, V), we may write 



X(po,p) = X 1 (p ,|p|)Y 1A1 (O p ), *(po,p)= XiCpo.lpDY^fip). (33) 

L=0,l,2 



The corresponding equations for the radial functions can be found by a partial wave decomposition of the kernel in 
eq. ( pS) ) and by carrying out the angular integration. An example of the system of coupled equation for the radial 
amplitudes in the case of one scalar exchanged boson is given in the ref. . In what follows the notation for the radial 
amplitudes are kept as in eq. (p2|) with the lower index indicating the value of the angular momentum L in eq. (p3|) . 



B. The transformation properties of the partial amplitudes 

Due to the parity invariance, (|3l|), only eight radial components are relevant to describe the deuteron amplitude, 
namely 

Pi, A?, A , A 2 , Vi, T°, To, T 2 . (34) 

Analyzing the behavior of the amplitude under the symmetry transformations, one can establish the properties of the 
components ([34]) . The invariance of the BS equation under the time- reversal operation T 

?Xm (PO, P) = 7 1 7 3 Xm * (PO, -Pb 1 7 3 (35) 
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and the complex conjugation JC 



iCx^(p) = (-i) m x d m (p) 



(36) 



imply that the seven partial amplitudes Pi, Vi, Aj, Ao,2, To,2 are real functions while the amplitude is purely 
imaginary, i.e., \ ( — — I \ . 

The Pauli principle implies that the amplitude xd(p) changes the sign if two nucleons are interchanged, i.e., 



Xd(po,p) = -Xd(-Po,-p)- 



(37) 



^From eqs. (37) and ( pq ) follows that A? and are odd functions with respect to the operation H(po — > — Po) 

nA?(p ,p) = -A?(pQ,p) ! nT°(p ,p)=-T?(po,p) (38) 

and the remaining six amplitudes are even functions of pq. This symmetry property is useful for the classification of 
the amplitude according to two-nucleon states with a given relative energy, i.e., the p spin classification. 

Table |l] summarizes the properties of the partial BS amplitudes in the representation (32) under the symmetry 
transformations. 



C. Observables 

Relying on the symmetry properties of the partial amplitudes, defined by eq. ([52]), the BS equation (p8|) has been 
solved numerically jj| for the deuteron at rest by performing a Wick rotation (po — > ipi)- In our present calculations 
we include six meson exchanges, 7r, to, p, er, rj and S, which describe the effective NN-forces. The set of the meson 
par ameters, such as masses, coupling constants and cut-off form factors, has been taken essentially the same as in ref. 
|L4j , where it has been obtained from a fit of the phase-shift analysis of the N N scattering and the binding energy of 
the deuteron. 

The BS amplitude does not have a direct probabilistic interpretation as the Schrodinger wave function. Moreover, 
there is no simple way to compare these two objects describing the same system, the deuteron. In order to make a 
comparison possible, we can compute the same matrix elements of observables in two approaches and compare these 
observables. 

For example, the normalization condition ( |30| ) in the rest frame of the deuteron and for p = 0, (D\N {0)7° N (0)\D) = 
2Md, is simply a charge of the deuteron associated with vector current. In the Wick rotated system (p — * ipi) and 
in terms of the partial amplitudes (BJ) it reads 



Md = 2 J d 2^L{-M D (P? + Af + 4Tf + V?; 

+ (2m N - M D ) (x o + 2 + - (2m N + M D ) [xf + xf^j 

+ 2 ~^Pi (X + - V2X+ + X - - v/2X 2 -) 
- (V2X+ + X+ - V2X - - X 2 ") }, 

where X* = y2(T ± A/2). Now we define the charge density p c h(|p|) as 



(39) 



^- {D \N(0h°N(0)\D) = J ^ p,(p4, |P|), (40) 

OO 

Pch(| P |)= J ^Pch(P4,|p|). (41) 

— 00 

This already may be compared with the corresponding nonrelativistic analogue, i.e., the square of the deuteron wave 
function in the momentum space, oc u 2 {p) + w 2 (p). 

In the same manner also the nucleon spin-density may be defined as density of the axial charge 



G 



^i-(^|iV(0) 757 iV(0)|^) = I dPi ^ 1 Pspi «(p4, |p|), (42) 

OO 

Pspin(|p|) = y^W„(p4,|p|)- ( 43 ) 

— OO 

In the nonrclativistic limit this density reflects the contribution of the D-wave admixture in the deuteron, oc u 2 (p) — 

Results of numerical calculation of the denned densities together with a comparison with their nonrelativistic coun- 
terparts obtained within the Bonn and Paris potentials are presented in figs. |lf and [|. All curves exhibit qualitatively 
similar shapes and are identical in the nonrclativistic region |p| < 0.5 GeV/c. If the momentum |p| increases, the 
deviations of the rclativistic results from the nonrelativistic ones becomes more significant, but still too small to be 
attributed to the relativistic effect. Rather it is compatible with the model differences. Particular attention is to be 
paid to the fig. ||, where the spin density is depicted. This function is rather sensitive to the internal spin-orbital 
structure of the deuteron. The fact that the "elementary oscillations" of the spin density in the potential models are 
reproduced by the solution of the BS equation might be interpreted as the relativistic structure of the deuteron is 
governed by the nucleon interaction in states with a positive energy and L = 0, 2, i.e., by 3 Si and 3 D\ configurations. 
Therefore, in spite of the quadratic forms of the partial amplitudes is not diagonal in eqs. ([59|) and (f43f), one can define 
the relativistic analog of the probability of the D-wave admixture in the deuteron. Carrying out the \p\ integration in 
eq. d43| ) and equating the result to (1 — 3/2 Pu) we find Pd ~ 5% (cf. [[l4|]), which is compatible with the probabilities 
of the Bonn (P D = 4.3 [§) and Paris ( p D = 5.9 @) potential models. 



D. The BS amplitude in different representations 



To have a closer analogue with the nonrelativistic consideration it is convenient to use another basis set of matrices 



in the decomposition of the BS amplitude. In the literature the two-spinor basis |25| is frequently used, which means 
an outer product of two spinors, representing solutions of the free Dirac equation with positive and negative energies. 
This basis is labeled by the relative momentum p, the helicities A^ and the energy spin pi of the particles p4[ , 
sometimes also called (J, Ai, A2, pi, P2) representation. In this case one usually adopts for the partial amplitudes the 
spectroscopic notation 2S+1 L^ 1,P2 , i.e., 

*s++, 3 sr, 3 Di + , 3 d^-, 1 Pt~, 'pr + , 3 ^+-, 3 ^r + . (44) 

Sometimes it is more convenient to change from the (J, Ai, X2, pi, P2) representation to the representation (J, L, S, p) 
where p is the projection of the total energy spin of the system. In this case the notation of the components is as 
follows 

Y T = {v°,vlv e s ,v°,u + ,u-,w+,w-), (45) 

where u,v,w correspond to L = 0, 1,2 respectively and o or e mean the odd or even parity relative to the p-spin 
function; the lower indices s,t denote the singlet and triplet spin configurations respectively. According to eqs. (|37| ) 
and (|36l), the amplitudes v°, vf are odd and v% , v° are even functions of po- The partial amplitudes in the basis (|44|), 
( p5| ) are of a more familiar form and show a more transparent physical meaning since they may be compared with the 
deuteron states in the nonrelativistic limit. It is intuitively clear (see also figs. [l]and[2]) that the two nucleons in the 
deuteron are mainly in states with L — 0, 2 and with positive energy so that one may expect that the probability of 
states with negative energies and L = 1 in eqs. ( ^i[ ) - (^) is much smaller in comparison with the probability for the 
3 S^~ + and 3 D^ + (or u + and w + ) configurations. Moreover, it can be shown that the waves 3 S^~ + and 3 Df + directly 
correspond to the 5* and D waves in the deuteron, while those with the negative energy vanish in the nonrelativistic 
limit. 

The partial amplitudes (^J) are defined through the following decomposition of the BS amplitude 

Xc(P0,p) = X)^(P0'IpI) V m(p). ( 46 ) 

a 

where a — { J, L, S, pi, pi\ labels different states of the system; 4> a denotes the partial amplitudes in eq. (|45|), V^(p) 
are the spin-angular functions 
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V^(p)=* L £ (LmS S \JM)(\ Sl \s2\S S )Y Lm (i>)U P sl(jp)U p s 2 2 (-v) 



2 2 



(47) 



In eq. (^) the quantities Ug(p) are the free nucleon spinors; the explicit matrix form for the spin-angular functions 
VjU (p) i s gi ven m the Appendix I. 

In order to establish a connection between the representation (32) and the spinor basis we represent the 

Dirac matrices in cq. (55) as a direct product of Pauli matrices of the nucleon spin a and the p-spin 



XD (p) = ^0 [/P - 2i(<r-f „)] + ip 2 ® [IA° + (<x-V)] +p 3 ®( £r .A) + 2/<8>(<r-f) 



The last two terms in eq. (fig) may be rewritten as 



1 



(o-A) + 2J®(o-T) = _(/ + /9 J )®( C r.X + ) + ^(/-p d )®(<T-X-), 



72 



V2 V 



(48) 



(49) 



Then eqs. ( |48"1 ) and (f49|), together with the symmetry properties of our partial amplitudes listed in the Table |, show 
that the desired relation between the two representations appears as follows: 



X+ 3 D~ 



3 Q++ ^ , Y+ 3 c— 

3 Pf ~ T?, 3 P X ° - Vi, 1 P1 ~ Pi 1 P 1 ° ~ A° 



x- 



The relation between (|34|) and 



and Aj, Tj are related direct 



can be established exactly. The components being odd in the relative energy v°, 
y to each other via 



(50) 



whereas the remaining six components are connected via linear combinations. By representing these amplitudes as 
six-component vectors, Y T — v° 7 u + , u~, w + , w~) and ^f T = (Pi, Vi, Xq , X^~, X^, X^), the transition from Y to 
VP is provided by a unitary transformation Y = U^> (with det(U) = — 1, and U U T = 1) with the following explicit 
form of the transition matrix 



U= 



(51) 
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V 3 


2 


c 


i-\A+C 2 
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2v^ l-V!+£ 2 
3 C 


[2 


2 


W1+C 2 


1+Vi+C 2 


2^2 1-\A + C 2 





V 3 


v's 


3C 
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3 C 
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[2 

V 3 





2v / 2 1 -V 1 +C 2 
3 C 


i+Vi+c 2 
c 


i-\A+C 2 
3£ 


2 

n/3 


[2 

V 3 


2^2 l-Vi+C 2 
3 C 





i-Vi+c 2 

3C 


i+\/i+C 2 
C 



V 

with £ = |p|/m. In the nonrelativistic limit, where ^ <gc 1 , the matrix U becomes diagonal 

J7 = diag(-1, 1,1, 1,-1,-1), 



(52) 



and our representation coincides with the one in the spinor basis. In what follows all formulae will be derived in terms 
of the partial amplitudes (H) or ( p5[ ) , nevertheless the numerical calculations are performed with our solutions (|34|) 
by utilizing eqs. (5^) and (pl|). 

Coming back to the normalization condition it is easy to show that eq. ( |39| ) may be transformed to a diagonal form 

•dp 4 d|p||p| 2 



2 



(2tt)4 



-(Y + (p 4 ,|p|),«y(p4,|p|)) = i, 



(53) 



which is exactly the normalization condition used in ref. JlJ] . In eq. (^3|) Y denotes the eight component vector (|4^) , 
and Cj is a diagonal matrix 
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uj = -diag(Af D , M D , M D , M D , M D - 2E p , 2E p + M D , M D - 2E p , 2E p + M D ) , (54) 

so that the integrand in eq. ( |53| ) consists of a sum of quadratic terms of radial functions Y a weighted with u a . Therefore 
each term, after integration, may be interpreted as pseudo-probability of finding the corresponding relativistic state 
in the deuteron. The result of our numerical calculations of the pseudo-probabilities is presented in Table ||. It is 
seen that an admixture of the negative-energy amplitudes affect the contribution of the positive-energy states. The 
appearance of the negative contributions of waves with negative p-spin is not a surprise; it follows from the physical 
meaning of the normalization condition according to that the contribution of each term in eq. (p3|) is the effective 
baryon charge in the corresponding state. The pseudoprobabilities of S and D waves (see Table 0) are close to the 
corresponding probabilities obtained in the nonrelativistic Bonn and Paris potentials, which is to be expected, since 
the deuteron is essentially nonrelativistic system. 

To investigate the behavior of the partial amplitudes and their nonrelativistic limits, we employ once more the 
normalization integral (|30|), now however in the form of eq. (53). Then, similar to eqs. ( fid]) and ([43|), we define the 
following functions ip a depending upon |p| by 



#*(|p|) = ^2 Jd Pi Lu a \Y a ( Pi ,\p\)\yM D . (55) 

Thus ip a may be regarded as the absolute value of the relativistic wave function of the deuteron in the state a (for 
instance, a — 5 corresponds to 3 S^~ + configuration, a — 7 to 3 Df + etc., cf. eq. (^5|)). 

Figures (|) and (|J) display the behavior of the relativistic wave functions ipQ and ip2 (solid lines) versus the relative 
momentum p| in comparison with the nonrelativistic S and D waves. We conclude that with an accuracy of model 
ambiguities in the nonrelativistic calculations (cf., the difference between Paris and Bonn wave functions, i.e., the 
dashed lines in figs. (|J) and (Q)) the large relativistic components are close to their nonrelativistic analogues up to 
|p| ~ to. However, there is a distinctive difference in the shape of the D waves in the two approaches. Namely, the 
nonrelativistic functions change the sign in the region |p| ~ to, whereas the BS component does not do so (cf. the solid 
line labeled as BS-I in fig. ||). To understand this we tentatively introduce an auxiliary definition of the relativistic 
D wave which is just the difference between the integrand in the normalization condition and the contribution of 
the 3 S^ + component, i.e., we introduce in the definition of D the contribution of all the negative energy states: 

"02 ~ V w+ + w~ + u~ + .... In this case only two wave functions ipQ and tp2 determine the normalization of the 
BS amplitude, and the correspondence with the nonrelativistic limit becomes one to one. In fig. |4] the function ^ is 
labeled by BS-II, and it is seen that it displays a minimum in the same region as the nonrelativistic functions, i.e., 
it has the same shape as the nonrelativistic D wave. One observes that the nonrelativistic D wave already mimicks 
relativistic effects, so that in calculations of relativistic corrections to the nonrelativistic approaches an overestimate 
of the magnitude of such corrections may occur. For completeness, in fig. ^ we present the wave functions for L = 1; 
since the waves u~ and w~ are negligibly small, even in comparison with the waves L = 1, they are not presented 
here. 



E. The vertex functions 



In studying the nonrelativistic correspondence of the solutions of the BS equation it is convenient to work with the 
BS vertices G(p; P) defined as 



x{p;P) 



(pi + m) G{p; P) (p 2 + m) 
(p\ — m 2 )(p\ — to 2 ) 



(56) 



^From eqs. (|46| ) and ( |5q ) it is possible to find a partial decomposition for the vertex G(p;P). In doing so, one 
introduces the two four- vectors of on-mass-shell particles corresponding to the Dirac spinors in eq. (|47|), i.e., 



ki = (E p ,p), k 2 = (E p ,-p), E p = Vp 2 + m 2 , p=(po,p)- (57) 
Then in eq. (|5^) the inverse propagator of the nucleons may be represented in terms of the vectors by 

S-\l) = ^ +P ~ m = ^- [(k - m)SZ 1 (l) + (k + m)Sl 1 (l) 



S- 1 (2) = ^-p-m=^- [(k - to)5Z 1 (2) + (h + m)Sl 1 (2) 



(58) 







where 



(59) 



S±(l)= (^+p TE p ^ , 5±(2)= (Ml-po^E 



Because of 

S-\l)U^(p) = sign( Pl )S Pl (l)[/f (-p) 

S-\2)U^(-p) = sign(p 2 )5 P2 (2)(7f (p) (60) 
the partial decomposition of G(p; P) reads 

G M ( P ;P) = ^G Q (p ,|p|)V^(-p), (61) 

a 

hence the partial amplitudes and the vertex functions are interrelated via the following simple expression 

r>o,|p|)= S Pl (l)S P2 (2)G a (poA P \). (62) 



The relation eq. (|62|) implies that the BS amplitudes (|46|) have sharp maxima around po — 0, while the behavior of 
the partial vertices G a (p , |p|) is predicted to appear as smooth functions of the relative energy (see also ref. @). 

The behavior of the vertex functions is shown in figs. || and [7] for the configurations 3 S^ + and 3 D^ + as functions 
of the relative energy po and momentum |p| in the Wick-rotated system. It is seen that the dependence of the vertex 
functions upon the relative energy is weak, hence one may expect that the nonrelativistic and relativistic vertices at 
po = have similar structures as functions of |p|. ^From this observation and eq. ( |62|) we establish another relation 
between the BS amplitudes and nonrelativistic wave functions. Below, as an example, we show how one can obtain 
the relativistic wave function for the 3 S^ + configuration from the BS amplitude. The energy dependence of the 
component u + is factorized into two parts, namely a dependence on the scalar propagators ( |59"| ) and a vertex function. 
Then using the smoothness of the vertex as function of po we replace it by its value at po = multiplied by a smooth 
function of po, i.e., 

" (M ' |p|) = w-^-rfi = pF^na m 

with £(0, |p|) = 1, where the dimensionless function £(po, |p|) reflects the energy dependence of the vertex function. In 
view of the smooth behavior of the vertices as function on pp , one may replace this function by a constant, £(0, |p|) « £0 
(with £0 ~ !)■ Then in the normalization integral eq. ( p3| ) the integration over the relative energy may be carried out 
explicitly and the remaining part corresponds to the square of the nonrelativistic wave function, i.e., we define the 
nonrelativistic limit of the BS amplitude u + as follows 

^o(1p1)=^ + (0,1p]) ^^=^ . (64) 



Similar definitions, using eqs. (|5^), (B4j), (59ft and (p2|) are valid for other waves. The generalized relativistic S and 
D waves in this manner are displayed in figs. || and The actual calculations have been performed with £0 = 1. A 
comparison with the corresponding nonrelativistic wave functions at |p| — ► shows that, by choosing the parameter 
£0 = 1) we slightly (by about 10%) overestimate the relativistic functions (see figs. |6| and |^). It is worth stressing 
that in our solution of the BS equation the relativistic D wave does not change its sign in the interval up to |p| ~ 1.5 
GeV/c. This is the most essential difference between the relativistic and nonrelativistic approaches in this region. 
Therefore, it is expected that the relativistic corrections to physical quantities in the deuteron up to |p| ~ 1 GeV/c, 
are relatively small; to distinguish them one should either compute observables which are known experimentally with 
a very high precision and sensitive to the spin structure, or find special processes where the large components are 
suppressed and only the states with negative energies are relevant. 



IV. THE STATIC CHARACTERISTICS OF THE DEUTERON 



Let us calculate the static moments of the deuteron in the BS formalism. The conserved electromagnetic current 
of the deuteron (0) in terms of the BS amplitude is given by 
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(P , ,X'\J lt \P,X) = -ieN D J d 4 pTr{xy(p / ;P')T^q) X x(p;P)S F (p 2 r 1 } , (65) 

where Sf{p2) — P2 + m, p' = P + q/2, P' = P + q, No = l/(2ir) /2Md- The quantity L^ is the photon-nucleon 
electromagnetic vertex, which is assumed to be of the on-mass-shell form 

r„(g) = 7^1 (? 2 ) - ^w v fZtf), (66) 

where — 5(7^,7^], and Ff are isoscalar Dirac (Pauli) form factors of the nucleon with i*i(0) = Fz(0) = 1/2, 
k = p,p + fj, n — 1, and /i Pi „ are the proton and neutron anomalous magnetic moments in units of the nuclear magneton 
l/(2m). The gauge invariance of the electromagnetic current in the ladder approximation has been proven in ref. [T^ ] 
(see also @). 

Now we have to interrelate the expression for the static moments (|2^) and ( p7| ) , which are determined in the Breit 
frame, and the BS amplitudes, which are numerically obtained in the rest frame of the deuteron. This relation is 
given by the general transformation rules 

Xx( P ; P) = A(£)xx(£- l P; PcmJA" 1 ^), (67) 

X x(p'\ P') = hr\C)x\W\ P c .m.)A(£), (68) 

K-\C)S f {\p - Pr'HC) = S F {\Pc, n . - C-'p)- 1 , (69) 

where A is the operator for spin-^ particles corresponding to the Lorentz transformation P = £P cm ., P' = C~ l P c . m ., 

A(£) = M ° + P ^ (70) 
y/2M D (E + M D ) 



with the corresponding Lorentz transformation matrix C 

~Vv \ 

n n 1 

(71) 



C 



/V^TTj -jrj 
10 
1 

V -y/fj V^+v- 



The direction of the boost is supposed to be parallel to qz- Then after the Lorentz transformation of the integrand 
in eq. (|65|), the matrix element takes the form 

(p',a'|JmI^a) = 

-ieN D [ d 4 pTr(xA'b';P c . m .)f M ( g 2 )xA(p;P c ,n.)5 F (ip c . m . - p)" 1 ^- 1 ^)] 2 ) , (72) 



2 



where 



r,(q)=A(C)r^q)A(C) (73) 
and the variable p' is represented via p and q as 

p> = Cp 1 ^ = C(pW + l q) = C? v + l C q, (74) 

with components 

p ' = (1 + 2 V )p a - 2^^T+7 ] p z - M D r), (75) 

p x '=p x , p v '=p y , (76) 

p z ' = (l + 2 V )p z -2^j^T+^p° + M D ^j x /T+V ] . (77) 

Eq. ( [72| ) is the starting point in evaluating the static moments of the deuteron in the BS formalism. The main 
peculiarities of this matrix element, in comparison with the familiar nonrelativistic expression, come from the Lorentz 
transformation and from the relativistic nature of the BS amplitude itself and might be summarized by 
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(i) effects of the negative-energy partial states (especially nondiagonal expectation of the current between 3 Si + 
and ipW'W 3pW.Ce) partial stateg)) 

(ii) a dependence of the amplitude upon the relative energy po ^ 0; in studying the static characteristics of the 
deuteron this effect is called retardation in the BS amplitude, 

(iii) an effect of boosting on the internal space-time variable, that is, the effect of C 7^ 1 

(iv) effects of the deformation of the BS amplitude concerning the booster A(£) 7^ 1. 

In fact, in the matrix element (|72|) these boost effects reduce to a deformation of the photon-nucleon vertex eq. (J73 

r a — 1 c rw* Tn , _^ n (cno ami doT^ onrl ifFftW W fV,o one ("Foh c^A ^7% i+ mo 



and to corrections from [A . In our case, i.e., as rj — * (see eqs. ( |2q ) and (p7|)) for the eqs. ( |72| ) and (73) it may 
be written: 

[A-W^ 1 + ^7073 + \ (78) 

A(£) 7o A(£) = 70 (79) 

A(£) 7l A(£)= 7l [A(/:)] 2 (80) 

A(£) la qA(£) = la q, (a = 0,1) (81) 

In what follows, the deviation of the quantity [A _1 (£)] 2 from unity in the matrix element eq. ( |72| ) we call the effects 
of the Lorentz boost in the BS amplitude. 



A. The quadrupole moment 



1. General formulae 

Accordingly to the eqs. (|27|), (|6q) and (|7|)-(^l]) the result for quadrupole momentum is presented as follows 

a,a' p,p' 

= EE [(*' P '\Qc\a p ) + (a /P '\Qh B \a") + (a'"' ' \Q M \a») + (a' p ' \Q L M B \a")] , (82) 

a, a' p : p F 

where the subscripts C and M mean the corresponding contribution of the charge and magnetic part of the photon- 
nucleon vertex (J6q) , and the superscript LB is the contribution of the Lorentz boost ([A(£) _1 ] 2 — 1). 

The corresponding matrix elements of the zero component of the deuteron electromagnetic current in the definition 
© take the form (A = 0, 1) 

J^ X) {P,P) = ^J -^ I Tr|xA(p';F)70XAb;P)^(|- J3 )- 1 |, (83) 
4 X - X) (P,P) = J ^Tr{xA(p';P)70XA(p;P)M^ -PrVft}, (84) 



e k 
2Mn 4m N 



Tr XA(p';P)A(£)( 7 og-<7 7 o)A(£)xA(p;P)^(--p)- 1 , (85) 



i(27r) 4 t ^-v~mi«m , w -^wr/"^ 2 

J^ A) (P,P) = -^^x (86) 



7 z(2tt) 4 Tr |^ A ^' ; P )^ 0< ? ~ <?7o)xa(p; - p)~ 7o73 j 



(87) 



As next step the partial wave decomposition of eqs. (p3| ) - ( p7| ) has to be performed. Then one expands the integrands 
in Taylor series around 77 = and carries out the limit 77 — > 0. It is clear that one has to keep corrections including 
0(77) in the wave function x(p ,p';P) and the matrix A(£). 

This scheme of calculation allows to investigate separately the contribution of different relativistic effects mentioned 
above. The eqs. Q) and (]87|) are new contributions which account for the effect of the boosted photon-nucleon vertex. 
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Moreover, the Lorentz deformation effect of the BS amplitude also is taken into account in these matrix elements 
through the relative momentum p' . 

Obviously, the main contributions to the quadrupole moment comes from the charge part (a' p |Qc|& p )) computed 
with the large S and D components of the BS amplitude. For these states, with p = p' = +1, one can recover the 
nonrelativistic formula for the quadrupole moment of the deuteron and separate the corrections due to the relativistic 
Fermi motion of the nucleons and the retardation in the relative energy 



Qd= E ( a ' + \Qc\a + ) = Q { p 



(+,+) i 

1 ^PO > 



(88) 



where the two terms in r.h.s. of the eq. ( pq ) reflect the existence of derivatives in respect to the momentum |p| and 
the relative energy in the corresponding integrands: 



and 



2M 



D 



-co 



dp \p\ 2 d\p\ M D . 



(89) 



1 {E^-mf, 



12 |p| 2 £ 2 



1 14E 4 , + bElm 2 - 3m 4 + 20Egm , „, 

—w(Po,\p\) 



120 



IpI 2 ^ 



+ Y^w + {p , |p|) 



V\W\ w+{po ' |p|) + ^ w+bo ' |p|) ^ +(po ' |p| 



y/2 3m 4 - 4E 4 , + 5££ m 2 + 5E 3 m 



60 

2£ p + 3m 



IpI 2 ^ 



"" + (Po, \p\)w + (p , |p|) 



i a 



p "-u + (p ,\p\)-^-^—w + (p ,\p\) 



20 E 



\/2 2E p - 3m 
"20 El 



|p|9|p| 
1 d 

w + {p Q , |p|) p ^| u+ bo, |p|) 



\/2 9 2 \/2 

+— | p |)_ tt +(p Dl |p|) + — U + ( P0 , IpI)^- 



r w + (p , |p|) 



5M D 



3 1 . ) 1 3^/2 

on2 W+ ^ ' l p l) 2 + |P|) |— r q]-, ™ + (P0, |P|) + y-rjU+ipo, \p\)w + (p , |p|) 

2 |pr IpI o\p\ \p\ 



+V2w + ( Po , \ p \)-LJL u + (po, |p|) + V2 U +(p , |p|) J- (p,,, | P |) 



IpI^IpI 



IpI^IpI 



+00 +00 



-V2 



2M 



D 



-co 
d 2 



d Po \p\ 2 d\p\ 1 | P | 2 m d 



<9 2 



(90) 



u+ bo, | P |)^-3w + (po, |p|) + w + {p Ql \p\)—^u + (po, |p|) 
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w + bo, |p|)^-3w + (p , |p|) 

+00 +00 

dp |p| 2 cf|p| 3 



2M T 



2 P0 



Md 



-00 



i{2ir) 4 10M D K M, 



D 
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V2 



1 + ir ) u+ (po> \p\)w-™ + (po, IpI) + ( i + ) w+bo, \p\)^-u+( POi |p|) 



9p 



dpo 



+w + (po, \p\)^—w + (po, |p|) 



+00 +00 

f f dp |p| 2 c?|p| |p| 



2M t 



i(27r) 4 5M D 



ci 2 P a \tw Md , ^ 
(1-^)(^p- — +Po) 



v/2 



<9 2 <9 2 

m + (po, IpI) a si i ^ + (po, IpI) + w + (p , |p|) Q M , u + (p , |p|) 



9p 9|p| 



9p 9|p| 



+w + (p , |p|) 



c) 2 



9p 9|p 



■w + (po, |p|) f , 



where u + (po, |p|) and u; + (po, |p|) represent the radial function of the corresponding partial state 3 S^~ + and 3 D* + . In 
the nonrelativistic approximation, E p — ► m^, pq/Mjj — ► 0, eq. j89| ) yields 



)(+,+) 



+00 +00 



2M D 10y/2 
8 2 



00 
5 d 



d p \p\ 2 d\p\ 

i(27r) 4 



d 2 



1 9 



a| P | 2 | P d|p| 



u + (j>o, IpI) 



w + (p , IpI 



3 



d| P | 2 |pd| P | |p| 

+00 +00 



^ + (P0,|P|) 



+ 



1 



2M B 20 



-00 



U 2 



d p Q \p\ 2 d\p\ 

i(2n) 4 [\d\p 



2 5 6 

|pqp| IpI 2 



X«,+ (po,|p|)(^-^). 



(91) 



The expression eq. ( Bll) h as not yet a "true" nonrelativistic form because of the integration over p . However, by 
making use of the eqs' (p3|), (Q) with £0 = 1 an d carrying out the po-direction integration explicitly, the familiar 
nonrelativistic expression |26| for the quadrupole moment is reproduced exactly 



ID 



1_ f d\p\ 
20 J (2tt)S 



|p|2 #qM#2([pD + gipi^dpD.^dPl) 



d|p| d|p| 



d|p| 



(92) 



IP|2 (*m) 2 +6MpI)) a 



where -0o(|p|) and ip2 (|p|) are defined by eq. (J64j) and correspond to the nonrelativistic S and D components of the 
deuteron wave function (see, also figs. || and UJTAs seen from eq. ( |92] ) the main contribution to the matrix element 
( p9| ) is expected to come from the interference of the positive S and D-states in the deuteron; the remaining terms 
with negative p- spins are the contribution of the relativistic Fermi motion. 

The second term Q^' + ^ in eq. (|88| ) and the matrix element of the Lorentz boost operator (|82] ) are of a pure 
relativistic nature and reflect the relativistic corrections to the quadrupole moment. For instance, for the positive 
states the corrections Q^ B ', <3^ +) = ( a ' + \Qc B \ a+ ) are 



a,a'=S,D 
+00 +00 



)(++) 
''LB 



2Mr 



dp \p\ 2 d\p\ ,v M D 2p 1 1 



i(2n) 4 



M D 'bM D E p 



(93) 



6E 2 , - 2mE p - m 2 



2 w+ (po, IpI) 2 + V2u+( Po , \p\)w+(p Q , |p|) 
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V2|p| 



d d 
u + (Po, IPDgipT w (Po> IpD + w+ (Po, l p Daj^| u+ ( po ' l p l 



d 



+ |p|w + (po, IPDgjpT w+ (Po, IpD 



+00 +OO 



dpo|p| a d|pL_ M D , ,1 |p 



2\f L) j J z(2^) 4 ( Ep 2 +Po hMr,E t 



(94) 



V2 



+ w + (po, \p\)—w + (p , |p| 
cpo 

After integration by part in Eq. ( p3| ) and (|94|) one obtains 

+00 +00 



d d 
u + {Po, \p\)jt- w + (Po, |p|) + w + (p , \p\)t^—u + ( Po , |p|) 



2M D J J i(2n 
-00 



d Po \p\ 2 d\p\ ^ M D 2 P0 2 1 
i — ( S p ^~ +Po)(l ~ x? - )?^— x 



(95) 



2^ _ m E p - m 2 



V2u+{p 07 \p\)w+{p , |p|) + \w+{p ai \p\f 



+OO +OO 



2Mr 



-00 



d P0 \p\ 2 d \p\ 1 |p| 

i(27r) 4 



5MlE p 



E P J 



(96) 



^ U +(po, |p|)w + bo, IpD + 2 w+ (po, IpD 2 



and vanishes in the nonrelativistic limit. 



It is seen that the magnitude of this term is of order ~ (jgL-/v;> 

In order to achieve self-consistency in the nonrelativistic approach to the deuteron form factors and electrodesinte- 
gration reactions, various relativistic corrections to the matrix elements must be taken into account, such as meson 
exchange currents and pair term contributions JT^j2^^8|| . In the covariant description of the deuteron these effects 
are partially accounted by calculating transitions between states with negative energies; the contribution of P states 
in the deuteron electromagnetic current corresponds to diagrams with nucleon-antinucleon pair creation in the old 
fashioned perturbation theory. Moreover, in ref. |l7j it has been shown that, considering the deuteron electrodesin- 
tegration process within the light-front dynamics, beside the dominant contribution of expectations with S and D 
waves, an extra matrix element with transitions between positive and negative energy states is relevant to describe 
the electrodesintegration amplitude. It has been also shown that the contribution of this extra component exactly 
reproduces the pair term corrections in the nonrelativistic limit. An investigation of the correspondence between the 
light-front dynamics approach and the BS amplitude has shown jl9| that the extra component in ref. |l7j may be 
imitated by transitions between a linear combination of the P-waves and S 01 D waves. Hence in our calculation the 
pair terms are taken into account via calculations of off diagonal expectation values of the relevant current between 
the S and P partial states (see also discussions in refs. ^,0). A more detail analysis of the nonrelativistic limit of 
the expression for the quadrupole moment with keeping leading corrections ~ 1/m will presented elsewhere. 



2. Numerical results 



The total expression for the quadrupole moment consists of a multitude of terms likewise eq. (89) with quadratic 
combinations of partial states and terms with second derivatives 9 2 /9|p| 2 , d 2 jdp\ and mixed ones d 2 /d\p\dpo com- 
puted between different partial BS amplitudes. Their analytical form has been evaluated by an algebraic formula 
manipulation code. Concrete numerical calculations have been performed by using our solutions of the BS equation 
for the partial amplitudes eq. (|34|) and the relations (|50|), (51). We find that the main contribution to the deuteron 
quadrupole moment gives the first term in eq. 



3) and that the transitions between energy even-states dominate, i.e., 
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Q++ = 0.2690 fm 2 . (97) 

This contribution is below the experimental data Qd — (0.2859 ± 0.0003) fm 2 [^9| by about 6%, nevertheless it is 
larger than the usual nonrelativistic calculations. This is an understandable effect because of the specific feature of 
the solution of the BS equation for which the sum of the pseudo-probabilities of the positive S and D waves is larger 
than 1. In this context, since the pseudoprobabilities of the remaining configurations are negative, the transitions with 
P waves are expected to play an important role in studying the static characteristics of the deuteron. A particular 
interests present calculations of off-diagonal expectation between the S and P partial states, which is predicted to 
replace the meson exchange contribution in nonrelativistic calculations ]l7| , ^9[ . Indeed, our numerical result points to 
a significant contribution of the mentioned matrix elements in comparison with other nondiagonal transitions 

(u + \Qc\v e s ) = 0.0052 fm 2 ; (u + \Qc\v°) = -0.0027 fm 2 (98) 

(for example, among other nondiagonal matrix elements the largest one is (w + \Qc\v 1 ;) = —0.00007 fm 2 ). 

The part of the quadrupole moment with odd (diagonal and nondiagonal expectations) gives a small negative 
contribution to the first term in eq. (^): (\Qc\)odd — — 0.0007 fm 2 Gathering together all the contributions we obtain 
Qc = 0.2706 fm 2 . 

An estimate of the corrections owing to the dependence on the relative energy, eq. J8S|), shows that they are rather 
small: Q^+ +) = 0.0006 fm 2 . 

The Lorentz boost corrections have been calculated and they are found to be negative. Their total contribution is 
Q C B — — 0.0029 fm 2 which, together with Qc, gives the final result for the electric part of the quadrupole moment 
of the deuteron Q = 0.2683. An important moment should be stressed here. The contribution of the Lorentz boost 
terms with nondiagonal transitions between S and P waves are of the same order of magnitude as those in eq. (|98l ) 
but opposite in sign 

(v + \Q c B \v e s ) = -0.0053 fm 2 ; (v + \Q^ B \v°) = 0.0027 fm 2 . (99) 

Equations ( p8| ) and (|9^) show that the contribution of pair creation terms in nonrelativistic calculations is predicted 
to be negligibly small for the quadrupole moment and confirm the qualitative results obtained in ref. |lq| . 

Note that our classification of the matrix elements into the main part and Lorentz boost corrections (cf. eq. (|82) 
rather conventional and does not reflect directly the contribution of relativistic effects. However, by using eqs. 
@ and @ we may present our results in the form 

Qd = Qnr + SQrei. = (0.2690 - 0.0007) fm 2 , (100) 

where Qnr is determined by the large components of the BS amplitude and does not depend upon the derivatives in 
respect to the relative energy and upon the Lorentz boost effects; SQ re i. is the contribution of all the remaining terms 
and, obviously, is of a pure relativistic nature. It is seen that the relativistic corrections to the quadrupole moment 
are negative and reinforce the discrepancy, although their magnitude is rather small. A similar conclusion has been 
drawn in ref. JToj ] from a qualitative analysis of the deuteron moment within the BS formalism. 

Another source of the relativistic corrections is the contribution of the magnetic part of the effective current ( p2] ) 
which vanishes in the nonrelativistic limit. Our calculation show that its contribution to the quadrupole moment is 
negative too, ( \Qm \ ) = — 0.0005 fm 2 , so that our final result for the deuteron quadrupole moment is Qd = 0.2678 fm 2 , 
i.e., the discrepancy in Qd is about 6%. 



B. The magnetic moment 

1. General formulae 

Accordingly to the eqs. @, @ and (|7|)-(|l]), the result for magnetic moment can be written as follows 

[ID = M+ + Mi- + M2- + M3-, (101) 



where the matrix elements between states with positive energies in eq. ( 101 ) are labeled by the subscript + , and 
the subscript — means that the corresponding matrix element implements at least one wave with negative energy. 
These matrix elements reflect the relativistic corrections. In order to emphasize the nonrelativistic analogue of 
the magnetic moment in the expression for the /i + we subtract the corresponding nonrelativistic formula and the 
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remaining part we denote as the relativistic corrections due to the Fermi motion effects. Then the functions /i.. 
can be represented by 



3 1 

A*+ = (Mp + Vn){P u + + P w +) - + - 2^ Pw+ + R+7 *- 102 ^ 



Mi- - \(lh + Mn)(A| + Pv°) + \{Pv t + Pa*) + \{Pv t +P K ) + Rl-, (103) 

1 5 

("2- =-(jip + Hn)P u - + P u- + 2<^P + Vn)P-m- - J P w- + R 2-, (104) 

^- = E cM, < ( 105 ) 

a . 6 

where a = n + , u> + , n~ , w~ , 6 = uj, v°, vf, v°, and are the pseudo-probabilities of the corresponding partial state. 
In eqs. (102)-(104) the diagonal expectations between states with L = 0,2,1 are written explicitly; the off-diagonal 
contributions are included into the terms R and /i3_ , where 



1, 2m , „,+ m „ + m „ + , 2m, „ 
- A*p + M« - 1 + )#i i?2 - (1 )-P 

.-^ p + lln -l-—) Hl --H 2 --H 3 -4(1-^ 

v/2, m 
-_( Mp + Mn -l-- 



(^ + ^-1-^)^ , (106) 



1 2m 1 1 2m 






1 rr<>«? 

5 9 


2 e 

Z TT V t > V t 

5 10 


+ 2 - { H<>< + 2H<><)- 


l<" Vt + 




r . 4to 2 

+ V2(/ip+Mn 1+ M2 


)Hl tK 




+V2(/X P + Mr* - l)i?5 tVse 







(107) 



1 , 4m % .„- m .„- m 



3. 2m, „ v2 , m, „- ,„- 

+ 4(1 - ^)iV + ^(Mp + Mn - 1 + • (108) 

The quantities C a - h and iT? ' Q are given in the Appendix II. Now the nonrelativistic formula for the magnetic moment 
may be recovered exactly by rewriting the term /i + in the form 

fi+ = fi N R + A^i + , (109) 

where 

3 1 

= (Mp + Mn) - T^P + Mn - 2)-^ 

reproduces the nonrelativistic formula and the relativistic corrections due to the Fermi motion effects are 

A M+ =R+~ (ftp + /!„) (P u - + P w - + P vl + P v o + P v e + P v o). (110) 
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Finally, the total contributions to the deuteron magnetic moment read as 

M-D = Pnr + AjU, 
Afi = R + + A/i_ + ^ 3 _, 

A/X- = -(fl p + fin) [\{Pvl + Pv S ) + (Pv t + Pvo) + 2P U - + X -P w 

+ i(Pspe + P v o) + l(P lp e + Hp.) + P„_ - |P W - + P d _ 



i?2 



(111) 



(112) 



Numerical results 



Explicit numerical calculations give for the total deuteron magnetic moment the value [Lb — 0.856140 (e/2m) which 
differs from the experimentally known moment [i exv — 0.857406 ± 10 -6 (e/2m) jl6) by less then 0.15%. This result 
consists of the nonrelativistic contribution plus relativistic corrections listed below 

(i) the main correction to the nonrelativistic value of the magnetic moment fi^R = 0.850718 (e/2m) that comes 
from the transitions between positive energy states and P-states (uf , v°, vf , v°) (the term in eq. (Ill)); it gives 
//3_ = 6.099- 10 -3 and contains ~ 0.71% of the total magnetic moment; (ii) relativistic corrections from the expectation 
values of p ositi ve energy states of the Lorentz transformation of the intrinsic variables in the BS amplitude (the term 
P+ in eq. (Ill)), which is found to be negative, i.e., R + = —9.75 • 10~ 4 ; (iii) the term Ap_, and the sum of transitions 
between states with negat ive e nergy (u~,u>~), and transitions between P-statcs themselves, and part coming from 
normalization effects (eq. ( |110D ); this is a positive contribution with A/i_ = 2.99 • 10~ 4 to the total moment. 

An analysis of our numerical results obtained for the off-diagonal expectation values between the S and P partial 
states shows that in contrary with eqs. (p8|)-(|99|), the contribution of terms like pair-creation corrections in this case 
do not compensate each other and give a total contribution to the magnetic moment ~ 0.35%, which is almost 50% 
of the total relativistic correction. 



V. CONCLUDING REMARKS 



In this paper we have investigated in some detail the numerical solution of the Bethe-Salpeter equation |9J with a 
realistic one-boson-exchange interaction. Special attention has been paid to a study of the relation of the partial BS 
amplitudes to the nonrelativistic wave functions and to the covariant description of the static characteristics of the 
deuteron. In our analysis we consider various basises used in defining the partial BS amplitudes and the transition 
from one basis to another. The representation based on the complete set of the Dirac 7-matrices and their bilinear 
combinations is found to be extremely convenient in computing the deuteron observables and processes with the 
deuteron pC| since in this case the dependence on the kinematical variables is mostly included in the definition of the 
partial amplitudes (except for one spinor propagator, which usually appears when computing diagrams for concrete 
processes, see ref. 0) and the matrix structure of the corresponding matrix element is almost independent on the 
intrinsic deuteron variables. However, in this representation an analysis of the deuteron structure in terms of familiar 
S, D, etc. components and an investigation of the correspondence of the obtained results with their nonrelativistic 
analogues is straitened. For this sake it is more convenient to use the p spin classification of the amplitudes for 
which a physical treatment of results is easier. In order to combine the advantages of these two representations the 
corresponding unitary transformation has been presented explicitly, cf. eq. (|5l|). With this a hand, calculations of 
various processes can be performed easily in the basis of the Dirac matrixes and the final expression may be treated in 
terms of the p-spin partial amplitudes by utilizing eq. (|5l"|). This scheme of calculation has been employed in order to 
compute the pseudoprobabilities of different partial states and to find the nonrelativistic limit of the amplitudes. In 
Section III different methods of comparison of our amplitudes with the nonrelativistic S and D waves are presented. 
Apparently, the most appropriate way to define the nonrelativistic limits of the BS amplitudes is to use the relation 
(|64|), which is based on an analysis of the behavior of the BS vertex functions in dependence on po and |p| and 
on the nonrelativistic relation between the vertices and wave functions in the momentum space. Numerical results, 
displayed on figs. [8] and ^ show that the generalized BS wave functions ( |64| ) are close to the nonrelativistic ones only 
for moderate values of |p|, while a difference occurs for |p| > m. This means that for rough estimates of possible 
relativistic effects one may calculate the corresponding nonrelativistic expressions by utilizing the wave functions (pi) 
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instead of the nonrelativistic S and D waves. Obviously, for a consistent investigation of the relativistic corrections 
it is necessary to use the covariant calculations with complete BS amplitudes. 

We have investigated the quadrupole and magnetic moments of the deuteron within the BS formalism by computing 
in the Breit frame the matrix elements of the electromagnetic current of the deuteron. In our analysis we considered all 
the possible relativistic effects connected with the Lorentz transformation from the rest frame of the deuteron to the 
Breit frame and with the dependence of the amplitude on the relative energy pq. By utilizing results of the investigation 
of the properties of the BS amplitudes performed in the Section III and their nonrelativistic limits, the static moments 
of the deuteron have been presented as a sum of two terms, one of them having a direct nonrelativistic analogue, 
the other one being of a pure relativistic nature. We pay special attention to the contribution of the nondiagonal 
expectation values between S and P configurations which are thought to include into the relativistic calculations the 
effects of pair currents widely discussed in nonrelativistic theories. It has been shown that for the quadrupole moment 
the different partial transitions between S and P components possess a noticeable magnitude, however, their summed 
contribution is found to be negligibly small (see eqs. ( p8[ ) and (|9^)), whereas for the magnetic moment these matrix 
elements give almost 50% of the relativistic effects. We obtain a good description of the experimental data for the 
magnetic moment. As for the quadrupole moment the computed value is below the experimental data by about 6%. 
That indicates that even a consistent relativistic computation does not perfectly describe the data in the impulse 
approximation. Probably, an adjustment of the operator of the electromagnetic current of the deuteron is needed, 
e.g., by including additional terms not accounted for within the present approach such as meson exchange currents 
with two- meson exchange diagrams or A isobars |jic|] . 

VI. SUMMARY 

In summary an analysis of the properties of the partial Bethe-Salpeter amplitudes, obtained as numerical solution of 
the BS equation with a realistic interaction, has been performed. In order to compare relativistic amplitudes with the 
nonrelativistic wave functions a method, based on the comparative analysis of the observables, has been developed. 
The static characteristics of the deuteron, i.e., the quadrupole and magnetic moments, have been computed within 
the Bethe-Salpeter formalism with satisfactory accuracy. 
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APPENDIX I 

The matrix form of the spin-angular functions (p) , eq. ( [47| ) , may be obtained explicitly by replacing the outer 
product of the free nucleon spinors U£ (p) by their direct product, U[J* (p) <£> U^ T (— p). The BS amplitude takes then 
the form 

Xd(po,p)U c =X>«fr°>IPl) r A"f(P)^c, (H3) 

with r^(p) 

rfttP) = * L E (LmS S \JM)(^ Sl ^ S2 \S S )Y Lm (p)U^(p)U^ T (~p), (114) 

Sl S^ra 

where Uc is the charge conjugation matrix, Uc — «727o- 

One can exploit the p spin dependence and replace r^(p) = T"} plP2 (p), where 
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v/2£;p(m + 


EpY 


&2 — 





(115) 



with a € {i, S 1 , J}. 

The spin-angular structures for concrete partial waves are shown in Table III. Here £ M is the polarization vector 
of the deuteron with the components in the rest frame given by 



e+i = (-i,-*.o)/V2, e_i = (i,-*,o)/V2, Co = (o,o,i). 

and the four- vector — (®i£m)- 



(116) 



APPENDIX II 



Below we list the explicit form of the quantities G°' and H" ,a in eq. (101). By introducing new functions G 



the mentioned quantities are expressed as follows: 



T = U ,V 



r = u ± y s 



T = W ± ,V B S 



T = W . V° 



± e 
T = U ,V% 



T=U ± 1 V° 



± e 
T = W ,V t 



T = W , V? 



C T 



V6 to 
~12M 

V6 TO 
15~M 



G r = — — 



V3 TO 
12 M 



G T 



y/3 TO 
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6 - T ~JI 



G\ + AG T 2 - AG T 3 -Gl- AGl 
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12 ^ ^13 ^ ^14 + u 15 ' Lt 28 

±G{ 6 — G 17 + Gl — G T ± G27 
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r _ V6 
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=t( - T 23 ~~ ^24 "+" u 18 — ^19 ^ ^2! 



2' 1 



= TO 

3 M 



G25 + G*26 + Gl + -Gl + Gl 



\/3 „ T 

+ _ ^ _K(jr 22' 



-4- ^*nr 

± -^-«G 22 , 



where the G"' a are integrals of the form 



N J dp4|p| 2 d|p|i4i(P4,|p|)[5 4 r a (p 4) |p|)]ya'(P4,|p|), 

and Ai(p4, \p\) are scalar functions; i?i may be either a differential operator of the type d/dp4, d/d\p\ or a scalar 
function (p^ — —ipo), which are summarized in the following tabular form 
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Analogously, the functions H"' a (H"' a = Hf) are of the same structure as G"' a with 
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TABLE I. The deuteron partial amplitudes and their transformation properties. 
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TABLE II. The pseudo-probabilities of the partial waves in the deuteron 
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TABLE III. Spin-angular functions f %y for the deuteron channel. 
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Figure captions 



FIG. 1. The nucleon density in the deuteron computed within BS formalism in comparison with the nonrelativistic results. 

FIG. 2. The nucleon spin distribution in the deuteron computed within the BS formalism in comparison with the nonrela- 
tivistic results. 

FIG. 3. The momentum dependence of the 3 S± + component defined by eq. (^) (solid line) in comparison with the 
corresponding nonrelativistic wave functions with Bonn and Paris potentials (dotted and dashed lines, respectively. 

FIG. 4. The momentum dependence of the 3 Df + component. The solid line (BS-I) depicts the result of computation by (pB[); 
the dotted line (BS-II) includes the contribution of P- waves (see text); short- and long-dashed lines depict the nonrelativistic 
wave functions with Bonn and Paris potentials, respectively. 

FIG. 5. The momentum dependence of the P waves defined by eq. (p5h 



FIG. 6. The behavior of the vertex function G(po, |p|) for the 3 S^ + configuration in the deuteron in dependence on p4 and 

IpI- 

FIG. 7. The same as fig. [] but for the 3 D^ + configuration. 



FIG. 8. The nonrelativistic limit of the 3 S^ + component defined by eq. ( |64| ) (solid line) in comparison with the nonrelativistic 
wave functions with Bonn and Paris potentials (dotted and dashed lines, respectively. 

FIG. 9. The same as fig. M but for the 3 D* + components. 
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Fig. [I]. L.P. Kaptari, A. Umnikov.... Bethe-Salpeter Amplitudes. 
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Fig. ||. L.P. Kaptari, A. Umnikov.... Bethe-Salpeter Amplitudes. 
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Fig. ||. L.P. Kaptari, A. Umnikov.... Bethe-Salpeter Amplitudes. 
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Fig. H L.P. Kaptari, A. Umnikov.... Bethe-Salpeter Amplitudes. 
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Fig. 0. L.P. Kaptari, A. Umnikov.... Bethe-Salpeter Amplitudes. 
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Fig. ||. L.P. Kaptari, A. Umnikov.... Bethe-Salpeter Amplitudes... 
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